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ABSTRACT 

Far Ultraviolet Spectroscopic i5a:p/orer observations are presented for WD 0621— 376, a DA white dwarf 
star in the local interstellar medium (LISM) at a distance of about 78 pc. The data have a signal-to- 
noise ratio of ~20-40 per 20 kms~^ resolution element and cover the wavelength range 905-1187 A. 
LISM absorption is detected in the hnes of D I, C II, C II*, C III, N I, N II, N III, O I, Ar I, and 
Fe II. This sight line is partially ionized, with an ionized nitrogen fraction of > 0.23. We determine the 
ratio D/0 — (3.9 ± } o) x 10^^ (Sc)- Assuming a standard interstellar oxygen abundance, we derive 
D/H w 1.3 X 10^^. Using the value of A(H I) derived from EUVE data gives a similar D/H ratio. The 
D I/N I ratio is (3.3 ± 1%) x IQ'^ {2a). 
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possible D/H variations that may occur within the inter- 
stellar medium (ISM) of the Galaxy. 

We present the first measurement of deuterium ab- 
sorption toward the white dwarf WD 0621—376 (Z, b = 
245.41°, -21.43°; V = 12.0). The star is a hot (Toff « 
62,000 K), metal-rich DA white dwarf (Holberg et al. 
1993) with a smooth stellar continuum. It lies at a pho- 
tometric distance of ~ 78 pc (Holberg et al. 1998). In 
agreement with the lUE observations, no He II absorp- 
tion (e.g., AA992.3, 1084.9) is seen in the FUSE spectrum; 
therefore continuum placement for the observed D I lines 
in the FUSE band is not affected adversely by the far-UV 
He II series. 

This article is part of a series dedicated to the first 
FUSE measurements of deuterium along sight lines to- 
ward FeigellO (Friedman et al. 2001), WD 2211-495 
(Hebrard et al. 2001), HZ43A (Kruk et al. 2001), G191- 
B2B (Lemoine et al. 2001), BD -1-28° 4211 (Sonneborn et 
al. 2001), and WD 1634-573 (Wood et al. 2001). Moos et 
al. (2001) present a synthesis of these results. 

2. FUSE OBSERVATIONS AND DATA PROCESSING 

FUSE is a dedicated observatory providing access to 
the 905-1187 A wavelength range at high spectral reso- 
lution. The design and performance of the FUSE spec- 
trograph have been described by Moos et al. (2000) and 
Sahnow et al. (2000), respectively. Two sets of FUSE 
observations of WD 0621—376 were obtained in histogram 
accumulation mode: On December 6, 2000, 19 exposures 
totaling an exposure time of 8.3 ks were obtained through 
the large 30" x 30" aperture (LWRS), while on February 3, 
2001, 21 exposures totaling 9.8 ks were obtained with the 
medium 4" x 20" aperture (MDRS) . The data are archived 
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1. INTRODUCTION 

Deuterium measurements place fundamental constraints 
on Big Bang nucleosynthesis (BBN), the amount of 
baryons in the Universe, and the chemical evolution of 
galaxies (see Tytler et al. 2000). Deuterium is thought to 
be formed by primordial BBN, and has subsequently been 
processed in stellar interiors. Therefore, the abundance of 
D should decrease with time and local D /H measurements 
should provide lower limits to the primordial D abundance. 
By comparing local measurements with measurements ob- 
tained for distant intergalactic clouds toward quasars, a 
better picture of the chemical evolution of the galaxies can 
be drawn. A key step is to establish a robust set of mea- 
surements for environments within the Galaxy. This has 
not yet been done, and the picture that now exists is com- 
plicated by several recent studies. Linsky (1998) obtained 
a nearly constant ratio of D/H — (1.5 ± 0.1) x 10^^ in 
the Local Interstellar Cloud (LIC), while other investiga- 
tions (Laurent et al. 1979; York 1983; Vidal-Madjar et al. 
1998; Hebrard et al. 1999; Jenkins et al. 1999; Sonneborn 
et al. 2000) suggest variations amongst the different sight 
lines probing gas beyond the LIC. Possible causes of D/H 
variations are numerous (Lemoine et al. 1999), but none 
have yet been definitely implicated in causing the observed 
variations. 

The Far Ultraviolet Spectroscopic Explorer (FUSE) is 
capable of directly observing deuterium absorption in the 
Lyman series lines (except Lyman-a) at ultraviolet (UV) 
wavelengths shortward of 1200 A. Access to higher order 
Lyman lines allows a larger variety of environments (more 
distant, denser clouds) to be observed than is possible with 
the Hubble Space Telescope or IMAPS. Therefore, FUSE 
measurements are essential for a better description of the 
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in the Multi-Mission Archive at the Space Telescope Insti- 
tute under the observation identifications P1041501 and 
P1041502. To reduce the eff'ects of detector fixed-pattern 
noise, the data were acquired using focal plane split mo- 
tions, wherein subsequent exposures are placed at different 
locations on the detector. Unfortunately, at the time of 
the observations this mode was not fully operational and 
caused misalignments of different channels during some 
exposures. For the LWRS observations, 8 exposures were 
lost in the SiC 2 channels. For the MDRS observations, 
only the LiF 1 channels were properly aligned; for LiF2, 
9 exposures were lost (for 7 exposures the detector was 
off); for SiCl, 6 exposures were lost; and for SiC 2 only 
two exposures were aligned. Airglow emission lines were 
reduced along with the stellar spectra, but their intensities 
were minimized as most of exposures were obtained during 
orbital night (except for two exposures for each observa- 
tion). 

Standard processing with the current version of the cal- 
ibration pipeline software (calfuse vl.8.7) was used to 
extract and calibrate the spectra. This version of the 
pipeline does not yet correct for the astigmatism or flat- 
field, and therefore we required any absorption features 
to be present in at least two channels to be considered 
real. The extracted spectra associated with the separate 
exposures of a given observation were aligned by cross- 
correlating the positions of strong interstellar lines, co- 
added, and rebinned to a nominal spectral resolution of 
~15, 000 (« 20 kms~^). The interstellar lines were shifted 
to kms~^ in the laboratory rest frame (note, however, 
that lUE spectra gave a heliocentric velocity of 40.5 ± 0.5 
kms~^, Holberg et al. 1998). Including or removing the 
two daylight exposures in each observation did not change 
the strengths of the airglow lines, and therefore all the 
exposures were co-added to maximize the signal-to-noise 
levels (S/N). To maintain an optimal spectral resolution 
and information on the fixed-pattern noise, neither the in- 
dividual segments, nor the LWRS and MDRS observations 
were co-added together. 

We present reduced spectra for several different detector 
segments in Figure 1. The data reveal numerous stellar 
and interstellar absorption lines. Prominent interstellar 
lines are identified. 

3. ANALYSIS 

We used two different methods to derive the col- 
umn densities of the interstellar species observed toward 
WD 0621-376: (1) The apparent optical depth method 
(AOD, Savage and Sembach 1991) was used to derive 
column densities and to check for unresolved saturated 
structures within the observed profiles. A direct inte- 
gration of the apparent column density profiles, Na{v) = 
3.768 X 10"Ta(t;)/[/A(A)], over the velocity range yields 
the true total column density in cases where the absorp- 
tion is weak (r < 1) or the lines are resolved. Equiva- 
lent widths were directly measured by integrating the ob- 
served intensity. (2) A profile fitting (PF) method using 
the code Owens was also used to obtain the column den- 
sities of D I, N I, O I, and Fe II. The reader should refer 
to Hebrard et al. (2001) for a full description of the PF 
method, and we will concentrate on describing in details 
the AOD results. The two approaches are different and 



complementary. Some of the issues involved in using these 
two methods are discussed by Jenkins (1996). By us- 
ing both approaches, we have a consistency check on the 
derived column densities and systematic errors associated 
with the measurements of each particular method. All 
methods are subject to a systematic uncertainty due to 
the possible presence of a broad component in the instru- 
mental line spread function (LSF). Several tests (Wood et 
al. 2001; Hebrard et al. 2001; Kruk et al. 2001) show how- 
ever that the uncertainty due to this effect is negligible 
for weak metal lines of the type discussed in this paper. 
Consequently, for the profile fitting method, the LSF was 
modeled with a single gaussian component. Exemples of 
how the LSF handling in the profile fitting affects the col- 
umn density results are discussed in § 3.2. 

We adopted wavelengths and oscillator strengths from 
the Morton (2000, private communication) atomic data 
compilation. This compilation is similar to the Morton 
(1991) compilation with a few minor updates to the atomic 
parameters for lines of interest in this study. 

3.1. Apparent optical depth method 
3.1.1. Interstellar absorption 

The stellar continuum was simple enough near the C II, 
C II*, C III, N I, N II, N III, O I, Ar I, and Fe II features 
to be fitted with low-order (< 3) Legendre polynomials. 
Equivalent widths and apparent column densities are sum- 
marized in Table 1, where lower limits indicate that the 
line contains some unresolved, saturated absorption that 
cannot be reliably estimated with the existing data. Com- 
parisons of the data from multiple channels were used to 
check for fixed-pattern noise and systematic continuum 
placement problems. For the measurements the following 
segments were used: SiC IB (905 -992 A); SiC 2A (917- 
1005 A); SiClA (1003-1091 A); LiF lA (987-1082 A); 
LiF IB (1094-1187 A); LiF2A (1087-1182 A). The last 
column in Table I indicates which segments were used to 
study each line. In this column, the absence of a seg- 
ment indicates that there may be a detector defect at this 
wavelength (e.g., SiC 2 A has more bad pixels than SiC IB, 
especially in the range 960-965 A near the N I lines). Al- 
though the measurements were cross-checked and were in 
agreement in the different segments indicated in Table 1, 
we adopted the values of equivalent width and apparent 
column density from the segments having the best quality 
data. Exceptions are the SiC 2A LWRS and SiC IB LWRS 
and MDRS channels that have comparable quality data; 
the values for lines measured on these segments were av- 
eraged (all the measured values were within the estimated 
1(7 errors of the individual measurements). 

To check if the lines contain some unresolved saturated 
structure, the apparent column density profiles of two or 
more lines of a given species with different values of /A 
were directly compared (N I and O I). The N I transitions 
at 963.99 and 965.05 A are affected by an unknown source 
of noise, and this results in a column density that can- 
not be reconcilied with the values derived from other lines 
(see Table 1). This discrepancy is unlikely to arise from 
the adopted oscillator strengths, as studies of other sight 
lines do not show any particular problem for these transi- 
tions. Since we have information from different channels 
(see last column in Table 1), fixed-pattern noise also seems 
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unHkely. Contamination by unidentified photosphcric stel- 
lar lines is the most likely source of confusion since this star 
is metal rich. These two lines were removed from the anal- 
ysis. Of the remaining N I lines, only N I Al 134.98 suffers 
from weak saturation effects. Similar comparisons were 
made for O I, and the transitions at 948.686, 971.738 and 
1039.230 A contain unresolved saturated structures (see 
Table 1). We note that, while the total column densities 
of all the unsaturated lines of N I and O I agree remark- 
ably well within the la errors, the strongest N I AA953.655, 
1134.415 and O I AA936.630, 976.448 lines yield the low- 
est apparent column densities, indicating some weak sat- 
uration effects might be present. For the other species, 
assuming they follow a similar curve of growth, only C II 
and N II are likely to be influenced strongly by saturation 
effects. 

Both C III and N III are blended with their respective 

stellar features. Figure 2 shows a stellar model (Tcs = 
62,000 K; log.g = 7.0, logC/H = -6.0; logN/H = -6.6, 
preliminary results from Chayer et al., in preparation) that 
estimates the stellar contribution (~30% for C III, ~7% 
for N III). The stellar model was aligned to the observed 
stellar linos, and an offset of w -1-14 kms""'^ was found with 
respect to the interstellar lines. N III suffers from blending 
with interstellar Si II A989.87 absorption, which is clearly 
visible in Figure 2, but the combination of relatively low 
spectral resolution and the proximity of the N III stellar 
feature did not allow us to model this line. We examined 
the lUE spectrum of WD 0621-376 to check for Si II ab- 
sorption at longer wavelengths, but its poor quality did 
not shed any light on the amount of Si II along this sight 
line. Therefore, the derived equivalent width of N III must 
be considered an upper limit. 

The spectra also show the detection of a very weak fea- 
ture Fe II at 1144.938 A. Though its measurement remains 
uncertain, its identification is unambiguous as the radial 
velocity of this feature is in agreement with the velocities 
of the N I lines in the LiF IB and LiF 2A segments. 

The 3a upper limits for the equivalent widths in Table 1 
are defined as Wmin = 3aSX, where a is the inverse of con- 
tinuum S/N ratio and SX « 0.1 A, which is approximately 
the average FWHM obtained from the resolved features. 
The corresponding 3a upper limits on the column den- 
sity arc obtained from the corresponding equivalent width 
limits and the assumption of a linear curve of growth. 

3.1.2. Deuterium, absorption 

D I Lyman-/? is well detected in the FUSE spectra of 
WD 0621— 376 (see Figure 3). Figure 4 shows some evi- 
dence for tentative detections of D I A972 and D I A949. 
Polynomial fits to the continua are shown for different seg- 
ments along with their respective normalized profiles with 
optimized gaussian fits (see below). Except for segment 
SiC 2A, all the measurements for D I AA972 and 949 yield 
detections with less than 3a significance (see Table 2). 
For segment SiC 2 A, both the continuum placement and 
the line detection remain very uncertain and fixed-pattern 
noise may contaminate this feature. The measurements of 
these lines and the 3a upper limits on W\ and log N were 
derived following the same methods described in § 3.1.1. 

Figure 3 displays the D I Lyman-/? spectra in the differ- 
ent segments. We note that the H I airglow is relatively 



strong and broad (~100 kms^^) in the LWRS data. Its 
strength varies in the different segments; for these data, 
the airglow in segments LiF 2B and SiC lA affects the D I 
measurement. The MDRS observations are less affected 
by the airglow, but it is clearly present. The LiF2B D I 
data are more affected by the H I airglow than the LiF lA 
data. 

The measurements for D I are summarized in Table 2. 
The continuum placement is more complicated than for 

the species described in § 3.1.1, requiring higher order (3- 
5) polynomial fits to the continuum. The continuum fits 
in the LWRS LiF 2B and SiC lA segments are very uncer- 
tain due to the airglow contamination; the measurements, 
however, agree with the other measurements (Table 2). 
Generally good agreement is found between the MDRS 
LiF lA, LiF 2B and LWRS LiF lA data. Of these, the 
MDRS LiF lA is the least affected by airglow, so we adopt 
this result (log Af(D I) = 13.75 ±0.09 dex, la) for compar- 
ison with the other value derived below. 

The continuum was also determined using the stellar 
model described in § 3.1.1 for the MDRS LiF lA observa- 
tions: The top panel of Figure 5 shows the model over- 
plotted on the data. While Figure 5 shows a reasonable 
stellar model, there are still some problems as some por- 
tions of the spectrum are not so well reproduced by the 
stellar model. To estimate the contribution from the H I 
feature, the normalized profile (bottom panel of Figure 5) 
was fitted with 3 gaussian absorption components using 
an optimized routine (Howarth et al. 1996), where the 
FWHM and the centroids were free to vary; changing (in 
a reasonable manner) the initial guess of these parameters 
did not change the fit. This method does not attempt 
to produce any realistic physical quantities for H I, but 
just tests the effective continuum placement for D I, as 
optimized gaussian fits can usually reproduce complicated 
profiles (Lehner et al. 1999). The bottom panel of Figure 5 
also shows the same fit with 2 gaussian components to the 
H I feature (FWHMs and centroids were fixed to the values 
derived with the 3 components fit) to indicate the degree 
of contamination from H I to D I. This approach yields a 
logarithmic apparent column density of 13.84 ± 0.10 dex 
(la). Considering other segments gives similar results, but 
these are more uncertain as the airglow is stronger. 

The above estimates of A''(D I) agree within their la 
errors, though the difference in the absolute values shows 
some uncertainties caused by the continuum placement. 
We adopt logA^(D I) = 13.79 ± 0.14 dex {2a), obtained 
by doing a weighted mean of the apparent column densi- 
ties derived using the two different continuum placements. 
Table 3 summarizes the definitive results for the different 
ions (for N I and O I, the values are weighted means of the 
values listed in Table 1) 

As discussed above, the S/N levels for the higher D I 
Lyman series were not high enough to make precise esti- 
mates of A^(D I). However, their 3a upper limits do not 
siiggest that the D I Lyman-/3 line is strongly saturated, 
if at all. Comparing the apparent optical depths of D I 
Lyman-/? with the unsaturated N I and O I lines (see Fig- 
ure 6) indicates that the D I Lyman-/? optical depths for 
the two continuum placements are only as strong as the 
strongest N I and O I features, suggesting that the line 
might not be entirely resolved (see § 3.1.1). However, sat- 
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uration corrections should be small (< 0.05 dex according 
to the N I and O I lines) and are less than the 2a error 
bars on the measurement. 

3.2. Results of the profile fitting and comparison 

In addition to estimating column densities using the ap- 
parent optical depth method, two independent profile fit- 
ting analyses using voigt profiles and using the informa- 
tion simultaneously from all the windows were performed, 
leading to two sets of results, PF 1 and PF 2. In PF 1, the 
profile fitting code Owens was used to fit the data in 32 
spectral windows, including all available unsaturated lines 
from D I, N I, O I and Fe II (i.e., N I A1134.980, O I 
AA948.686, 971.738 and 1039.230 were excluded). In PF2, 
a total of 24 windows wore used including lines from D I, 
N I, and O I. For Lyman-/?, only the MDRS data were used 
in both cases because the LWRS data are significantly af- 
fected by airglow. The results for these two methods are 
listed in Table 3. Although there are differences both in 
the absolute values and in the estimated 2a errors, all re- 
sults agree within \a. Apart from the different number of 
windows used, the differences in the results arise from the 
strategies adopted to deal with the widths of the LSF. The 
LSFs arc not yet well characterized and are known to vary 
from one segment to another. In PF 1 the LSF widths of 
about 60% of the windows were fixed to a pre-defined value 
(ranging from 8.7 to 12 pixels, depending on the window), 
while in PF 2 the LSF widths were free to vary, with the 
constraint that they remain within plausible values (in the 
accepted fits the maximum FWHM was 15 pixels and the 
minimum was 6 pixels).^ 

The largest differences, found for the N I column densi- 
ties, arc duo to the fact that in PF 2 high column density 
values implying saturation of some of the fitted N I lines 
- were permitted by the adoption of broader LSFs. 

The top and middle panels of Figure 7 show the best 
fits to the D I Lyman-/? fine in the LiF lA MDRS chan- 
nel, as obtained by PFl (logiV(D I) = 13.89 dex) and 
PF2 (logA''(D I) = 13.86 dex), respectively. As expected, 
there is no difference in the fitted profiles, but the dif- 
ference in the adopted fit parameters can be seen in the 
profiles plotted before convolution with the LSF. The un- 
convolved profile is slightly narrower (and deeper) in PF 1 
than in PF 2 because of a slightly lower 6- value (7.1 vs 8.9 
kms""'^) required to compensate for a slightly higher LSF 
width (10.5 vs 9.4 pixels in the LiF lA MDRS data). 

The bottom panel of Figure 7 shows a fit performed with 
Owens using the PF2 method but fixing the D 1 column 
density to the best value (log A''(D I) = 13.79 dex) derived 
in § 3.1.2 using the AOD method. In this case, the un- 
convolved profile is shallower than the other two, but the 
convolved profile is not, because this is compensated by 
a narrow LSF (8 pixels). For the AOD value the optical 
depth of the Lyman-/? line is smaller than 1. For the two 
best fits in PFl and PF2, the Lyman-/? line is slightly 
saturated, explaining why the D I column density results 
obtained by the two profile fitting studies are higher than 
those obtained by the AOD method. 

The results obtained by the three independent investi- 

® Note that the dispersions of the FUSE spectrographs are ~6.7 and ' 

Meyer et al.'s value (O/H = (3.19 ± 0.14) x 10"*) was corrected for the recommended oscillator strength of O I A1356 (/ 
instead of / = 1.25 x IQ-^) by Welty et al. (1999). 



gallons are shown in Table 3. They all agree within the 
estimated la errors. Our adopted result, shown in the last 
column of Table 3, is the simple mean between these three 
results. 

4. NEUTRAL AND PARTIALLY IONIZED GAS 

O 1 is an excellent tracer of neutral gas as its ioniza- 
tion potential and charge exchange reactions with hydro- 
gen ensure that the ionization of H I and O I are strongly 
coupled. A^(0 I)/A^(H I) is therefore an excellent ap- 
proximation of the ratio O/H. (Similarly the ratio D/0 
is well approximated by A^(D 1)/N{0 I)). Therefore, us- 
ing the Meyer et al. (1998) interstellar O abundance'^, 
O/H = (3.43 ± 0.15) X lO-'^, and the O I column den- 
sity along the WD 0621— 376 sight line, we derive a to- 
tal H I column density of 18.73 dex. We note, however, 
that the Meyer et al. sight lines typically sample more 
distant (100-1000 pc) gas than the WD 0621-376 sight 
line. From the FUSE spectra, it is not possible to derive 
a reliable H I column density, as the H I transitions are 
essentially confined on the flat part of the curve of growth 
for this sight line. Previously, Holberg et al. (1998) quoted 
log A^(H I) = 18.15 ± 0.15 dex (la), which is not in agree- 
ment with our value derived above. Using EUVE spectra, 
Wolff et al. (1998) found 18.70 dex, in good agreement 
with the value derived from O I, but unfortunately no er- 
ror is quoted. Considering the lUE spectrum, we note 
that it is not possible to obtain H I from these data due 
its relatively low quality. Only a combination of a very 
high S/N spectrum and a very good stellar model would 
help to directly determine the H I column density to a 
reasonable degree of accuracy. 

With the derived N I column density and the Meyer 
et al. (1997) intersteUar N abundance, N/H = (7.5 ± 
0.4) X 10~^, one finds 18.47 dex for H I, a value of 
about 48% lower than the value derived from O I. The 
FUSE bandpass allows us to directly estimate the quan- 
tity of ionized nitrogen, as N II and N III are accessible. 
N III is too poorly determined to be of much use as the 
amount of Si II remains unknown. However, the ratio 
N I/(N I-l-N II) < 0.77 again suggests that a substantial 
fraction (>23%) of N is ionized. Therefore, the discrep- 
ancy in values of A^(H I) derived from N I and O I may be 
due in part to ionization of N in H I gas along the sight line. 
While the charge exchange with hydrogen should make N I 
a good indicator of the amount of H I along the sight line 
(Sofia and Jenkins 1998), Jenkins et al. (2000) noted that 
in regions where ^ nuo , N may show a deficiency of its 
neutral form, similar to what is observed for argon. 

A similar result holds for the ionization of argon. Using 
our derived Ar I column density and the Ar abundance 
from Anders and Grevesse (1989), Ar/H = (3.6 ± g;?) x 
10~^, we derive an H I column density of about 18.30 dex, 
which, when compared to the value derived from O I, im- 
plies that 65% of the argon is ionized (neither Ar nor N is 
depleted significantly onto dust grains). These results for 
both N I and Ar I indicate that photoionization may be an 
important process along this sight line. Towards four other 
nearby white dwarfs Jenkins et al. (2000), similar results 

^6.2 mA per pixel for LiF and SiC channels, respectively. 

1.16 X 10"^ 
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were found, favoring the idea that the LISM ionization is 
maintained by a strong extreme ultraviolet flux (yet to be 
identified) rather than an incomplete recombination of he- 
lium from a previous episode of collisional ionization at 
high temperatures. The presence of C III also implies that 
the gas is partially ionized. Fe II is largely deficient (—1.3 
dex), but this is at least partly explained by iron being 
depleted onto dust grains. 

5. d/o, d/n and d/h ratios 

Using the results from Table 3, and following the above 
discussion (D/O = D I/O I), we calculate D/O = (3.9 ± 
1\q) X 10~^ (2o'). Using the interstellar O abundance, the 
D/H ratio can be derived: D/H « 1.3 x 10~^. A similar 
D/H ratio (D/H - 1.4 x lO"'"^) is obtained using N(H I) 
derived from the EUVE spectra. 

The D I/N I ratio is (3.3 ± 1%) x 10'^ (2a). As there 
is a non negligible fraction of ionized nitrogen and some 
of this fraction can be either in H I or H II regions, it is 
not possible to derive the ratio D/H using the nitrogen 
abundance. 

Toward WD 0621— 376, the H I column density remains 
relatively uncertain, but D/H is consistent with the LIC 
ratio ((1.5 ±0.1) x 10^^, Linsky 1998) and also consistent 
with other measurements within the LISM (Moos et al. 
2001). Since in the FUSE bandpass, D I and O I can be 
generally measured simultaneously, D/O ratio (more pre- 
cisely measured) can be used instead of D/H. To be noted 
that if D/O becomes available in distant quasars, the D/O 
ratio would be very sensitive to astration since D I is de- 
stroyed while O I is enhanced through the various stellar 
cycles. Initial measurements of D/O ratios in the LISM 
(Moos et al. 2001) show no significant variations which also 



favors constant D/H and 0/H ratios in the LISM. Other 
results suggest variations at larger distances (e.g., Jenkins 
et al. 1999; Sonneborn et al. 2000), but more observations 
beyond the LISM wall are needed to have a better under- 
standing of these possible variations in the ISM. 

6. SUMMARY 

We have presented FUSE observations of the DA white- 
dwarf star WD 0621-376, located at about 78 pc. The 
spectra show interstellar absorption from D I, C II, C II*, 
C III, N I, N II, N III, O I, Ar I, and Fe II. Using both the 
N ISM abundance and the N I and N II column densities, 
we showed that N is ionized by a fraction >23% (Ar is 
ionized by about 65%). These results confirm the recent 
suggestion by Jenkins et al. (2000) that a diffuse ionizing 
flux (> 24.6 eV) creates the high fraction of ionized helium 
observed in the LISM. 

The D/O ratio along the WD 0621-376 sight line is 
(3.9 ± i q) x 10~^ {2a), implying a D/H ratio of approxi- 
mately 1.3 X 10~^. This D/H ratio is also supported by the 
H I measurement obtained from EUVE. The D I/N I ratio 
is (3.3 ± X 10-1 (2a). Both D/O and D/H ratios are in 
agreement with other LISM measurements, implying that 
D/O and D/H are relatively homogeneous in the LISM 

This work is based on data obtained for the Guaran- 
teed Time Team by the NASA-CNES-CSA FUSE mission 
operated by the Johns Hopkins University. Financial sup- 
port to U. S. participants has been provided by NASA 
contract NAS5-32985. French participants are supported 
by ONES. This work was partially done using the profile 
fitting procedure Owens. f developed by M. Lemoine and 
the FUSE French Team. 
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Table 1 

Interstellar absorption lines (but D I): Equivalent widths and AOD results 



Ions 


A 






r 




S/N*- 


(mA) 


(cm^^) 


Segment' 


CII 


1036.337 


1, 


.18 


A 


1 n~ 


-1 


33 


98.9 ±2.5 


> 14.19 


1,(2,3) 


CII* 


1037.018 


1, 


,25 


V 
A 


1 n~ 


-1 


54 


12.8 ± 1.5 


13.06 ±0.05 


1,(2,3) 


cm 


977.020 


7, 


.62 






-1 


24 


62.5 


> 13.14 


4,5 


NI 


953.415 


1, 


.32 


X 


1 n~ 

lU 


-2 


32 


21.5 ±3.0 


14.34 ±0.07 


4,6 




953.655 


2, 


.50 


A 


-LU 


-2 


30 


32.9 ±5.0 


14.29 ±0.07 


4,5,6 




953.970 


3, 


.48 


X 


1 n~ 
±u 


-2 


30 


41.2 


14.24 


4,5,6 




963.990 


1, 


.48 


X 


1 n~ 

iU 


-2 


37 


30.8 ±2.3 


14.47 ±0.04 


6,(4) 




964.626 


9, 


.43 


X 


1 n~ 

iU 


-3 


30 


16.4 ±2.9 


14.35 ±0.08 


5,6 




965.041 


4, 


.02 


X 


±u 


-3 


32 


16.0 ±3.3 


14.72 ±0.15 


5,6 




1134.165 


1, 


,52 


V 
A 


1 n~ 

-LU 


-2 


48 


34.1 ±3.1 


14.35 ±0.04 


8,(7) 




1134.415 


2, 


,97 






-2 


35 


51.6 ±2.8 


14.31 ±0.04 


8,(7) 




1134.980 


4, 


,35 


X 


lU 


-2 


44 


56.2 ±2.8 


> 14.18 


8,(7) 


Nil 


915.613 


1, 


,65 


A 


-LU 


-1 


20 


51.5 : 


> 13.74 : 


5,6 




1083.994 


1, 


,15 


X 




-1 


40 


58.2 ±2.8 


> 13.80 


3,(9) 


NIII 


989.799 


1, 


,23 


V 
A 


-LU 


-1 


21 


< 79.0 


< 13.94 


1,(2,3) 


01 


924.950 


1, 


,54 


X 




-3 


26 


17.8 ±3.9 


15.23 ± 0.08 


5,6,(4) 




929.517 


2, 


,29 


X 


10- 


-3 


25 


21.8 ±3.9 


15.22 ±0.08 


4,5,6 




936.630 


3, 


,65 


X 


10- 


-3 


26 


37.1 ±4.5 


15.20 ±0.10 


4,5,6 




948.686 


6, 


,31 


X 


10- 


'3 


23 


57.1 ±4.4 


> 15.13 


4,5,6 




950.885 


1, 


,58 


X 


10- 


-3 


21 


18.8 ±4.3 


15.22 ± 0.09 


4,(5,6) 




971.738 


1, 


,16 


X 


10- 


-2 


23 


69.6 ±4.7 


> 15.04 


4,5,6 




976.448 


3, 


,31 


X 


10- 


-3 


23 


34.2 ± 6.8 


15.17 ±0.15 


4,5,6 




988.773 


4, 


,65 


X 


10- 


-2 


25 


... / 


... / 


4,(5,6) 




1039.230 


9, 


,20 


X 


10- 


-3 


32 


79.6 ±3.4 


> 15.11 


1,(2,3) 


Si II 


989.873 


1, 


,87 


X 


10- 


-1 


26 


e 


e 


1,(2,3) 




1020.699 


1, 


,64 


X 


10- 


-2 


41 


< 7.3 


< 13.68 


1,(2,3) 


Ar I 


1048.220 


2, 


,63 


X 


10- 


-1 


41 


17.6 ±2.8 


12.86 ±0.07 


1,(2,3) 




1066.660 


6, 


,65 


X 


10- 


-2 


40 


... g 


... g 


1 


Fell 


1144.938 


1, 


,06 


X 


10- 


-1 


42 


7.3 ±2.5 


12.88 ±0.15 


7,8 



Note. — Uncertainties are Icr error. Upper limits indicate that no feature is 
present and arc 3cr estimates, except for N III. Lower limits indicate that the 
absorption line is saturated. Colons indicate that the value is uncertain. (See 
text for more details). 

(o) Rest frame vacuum wavelengths and oscillator strengths are from Morton (pri- 
vate communication, 2000). 

(6) 1: LiFlA MDRS: 2: LiF lA LWRS; 3: SiC lA LWRS; 4: SiC2A LWRS; 5: 
SiClB MDRS; 6: SiC IB LWRS; 7: LiF2A MDRS; 8: LiF IB MDRS, 9: SiC lA 
MDRS. Segments between brackets mean that the lines measured on those dif- 
ferent segments are consistent with the presented result, but was of lower quality 
and arc therefore not included in the final results. The S/N level was averaged 
when more than one segment was used (though generally these ratios were simi- 
lar in the different segments), (c) The stellar contribution (~30%; see text) was 
removed from these measurements. 

(d) N I A953.970 is blended with N I A954.104 making these values uncertain. 

(e) The stellar contribution (^7%; sec text) was removed from these measure- 
ments, but N III A989 is also blended with Si II A989.9. This measurement 
supposes that the line is not saturated. 

(/) This line is a triplet of lines blended at the FUSE resolution. 

(g) Ar I A1066.660 is blended with the strong stellar Si IV feature at 1066.6 A. 
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Table 2 

d i measurements: aod results 





/ 


S/N 






Segment 


A 






(mA) 


(cm-2) 




1025.443 


7.91 X 10-2 


35 


34.2 ±6.8 


13.75 ±0.09 


LiF lA MDRS 






35 


41.5 ±8.7* 


13.84 ±0.10* 


LiF lA MDRS 






27 


32.9 ±7.0 


13.72 ±0.10 


LiF2B MDRS 






44 


37.5 ±5.8 


13.77 ±0.08 


LiF lA LWRS 






26 


40.4 ± 12.7 


13.80 ±0.15 


LiF 2B LWRS 






23 


39.3 ±9.2 


13.81 ±0.16 


SiC lA LWRS 


972.272 


2.90 X 10-2 


20 


11.5 : (< 14.8) 


13.67 : (< 13.79) 


SiC IB LWRS 






20 


13.2 : (< 14.8) 


13.74 : (< 13.79) 


SiC IB MDRS 






23 


18.9 : 


13.92 : 


SiC 2A LWRS 


949.484 


1.39 X 10-2 


27 


10.8 : (< 11.3) 


13.98 : (< 14.00) 


SiC IB LWRS 


Note. — 


Uncertainties 


are 


la error. Upper 


limits in brackets 


arc 3(7 estimates. 



Colons indicate that the value is uncertain. (See text for more details). (*) Stellar 
continuum, see text for more details. 



Table 3 

Column densities with 2a errors 



Ions 


AOD 


PFl 


PF2 


Adopted 


D I 
C 11* 
NI 
01 
Ar I 
Fell 


13.79 ±0.14 
13.06 ±0.10 

14.33 ±0:04 

1 c o-i _|_ 0.10 

12.86 ±0.14 
12.88 ±0:26 


13.89 ± g:05 

14.31 ± 0.04 
15.29 ± 8:06 

13.00 ± 0:09 


13.86 ±0.09 

14.37 ±0.12 
15.27 ±0.09 


13.85 ±0.09 
13.06 ±0.10 
14.34 ±0:09 

15.26 ±0:08 

12.86 ±0.14 
12.94 ± 0:23 



Note. — AOD gives the results of the apparent optical depth 
method described in § 3.1. PF 1 and PF2 refer to the two indepen- 
dent profile fitting analyses described in § 3.2. 
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Fig. 1. — Far-ultraviolet spectra of WD 0621—376 with identification of interstellaj: lines. 
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Fig. 2. — Spectra of C III and N III (histograms) and a model of stellar absorption (dotted lines). The N III feature shows clearly that 
interstellar Si II A989.9 must also be present. 



Deuterium toward the white dwarf WD 0621—376 



9 




Fig. 3. — Spectra of D I Lyman-/3 in the different segments and local polynomial fits to the continua. Note how the shapes of the lines differ 
between segments in the LWRS data (even for essentially night data) because of the H I airglow, while for the MDRS the effect of airglow is 
smaller. 
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Fig. 4. — Detection (2— 3a) of the D I lines in the SiC channels. The left panel shows the calibrated flux in different segments (and 
wavelengths), while the right panel side presents their respective normalized profile with a gaussian fit. 
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Fig. 5. — Spectrum of D I Lyman-/? in the LiF lA MDRS segment. Top panel: calibrated flux with the stellar continuum. Bottom panel: 
the resulting normalized profile with optimized gaussian fits (the solid line is the 3 gaussian component fit (D I and H I), while the dotted 
line is the 2 gaussian component fit (only H I); see text for more details). 
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Fig. 7.— Spectrum of D I Lyman-/? in the LiF lA MDRS segment with fits from PF 1 and PF2 (see § 3.2). The fitted profiles as well as 
the continua arc drawn by a solid line while the dashed line shows the intrinsic profiles before convolution by the LSFs. The bottom panel 
shows the value derived with an Owens fit using the PF2 method but with the D I column density set to the value (13.79 dex) derived from 
the AOD method. The dotted lines shown in the three panels indicate the continuum shown in the top panel of Figure 3. 



